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The bounds on dark matter models have been updated with the first result of AMS-02 on positron 
fraction in cosmic rays. The channels of annihilation/decay into r+r~ pairs have been excluded 
by the Fermi-LAT Galactic diffuse emission data. The channel XX ~^ IJ'^fJ-~ is allowed only for 
very cored dark matter distribution model (e.g., isothermal-sphere profile) while x ~^ M^A*" i^ 
viable. Interestingly, at energies < fOO GeV a positron/electron component from either dark matter 
annihilation/decay or pulsar seems to be favored by the AMS-02 data. In the case of annihilation 
of dark matter with a rest mass ~ 100 GeV, the cross section into electron/positron pairs is 
10~^® cm'' s~^. We also suggest that the primary electron spectrum may display hardening at ~ 240 
GeV, as observed by ATIC, CREAM and PAMELA in the cosmic-ray-nuclei spectra. Depending on 
the degree of the spectrum hardening, the subsequent positron-to-electron ratio can get flattened 
or even decreased, which may account for the flattening of positron fraction at energies > 200 GeV 
discovered by AMS-02. 
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I. INTRODUCTION 

In the standard cosmic ray model, most cosmic ray 
electrons are from supernova remnants while the cos- 
mic ray positrons are mainly produced through hadronic 
processes as cosmic ray protons collide with intergalac- 
tic hydrogen [H, In conventional approach, the in- 
jection spectrum of electrons (positrons) is taken as a 
single power-law. Since diffusion and electron/positron 
cooling are more efficient in higher energies, the spec- 
trum should soften with energy and the positron-to- 
electron ratio (i.e., $e+/($e+ + ^e-)^ where $ is the 
flux) should drop with energy monotonously Q. Hence 
there should be no prominent feature at TeV energies in 
cosmic ray electron/positron total spectrum, neither in 
the positron-to-electron ratio. The situation has changed 
dramatically since 2008. The Advanced Thin Ionization 
Calorimeter (ATIC Polar Patrol Balloon-borne Elec- 
tron Telescope with Scintillating Fibers (PPB-BETS [l). 
The High Energy Stereoscopic System (HESS ^) and 
Fermi Large area Telescope (Fermi-LAT [7|) reported 
cosmic ray electron/positron total spectrum up to TeV 
and found a hardening/bump in the energy range from 
100 GeV to 1 TeV. Almost simultaneously the Payload 
for Antimatter Matter Exploration and Light-nuclei As- 
trophysics (PAMELA) discovered an unambitious rise 
of the positron-to-electron ratio above ^ 10 GeV Q. 
Such a peculiar rising behavior has been confirmed by 
Fermi-LAT though in an indirect way [9] . People call the 
well-established electron/positron total spectrum hard- 
ening/bump and the increasing positron-to-electron ra- 
tio, unexpected in the standard model, the excesses or 
anomalies. These interesting features draw a lot of at- 
tention, and various physical origins, in particular new 
astrophysical sources and new physics (dark matter), 
have been extensively explored [see [lol for recent re- 



views]. Very recently, the AMS-02 collaboration has re- 
leased their first result on positron fraction in cosmic 
rays, which confirms the positron excess with unprece- 
dented accuracy up to the energy ^ 350 GeV |11]. 

On the other hand, the protons and Helium are the 
most abundant cosmic-ray components, and the spec- 
trum of these cosmic-rays up to the so-called "knee" can 
be described by a single power law [l], Q . Surprisiiigly, 
the spectra of protons (Helium) measured by ATIC [l2] , 
Cosmic Ray^Energetics And Mass (CREAM [3) and 
PAMELA 14j show a remarkable hardening at the mag- 
netic rigidity ^ several hundred GV (GeV/nucleon). We 
call such a kind of spectrum hardening the cosmic ray 
nuclei excesses, which challenge the current paradigm of 
cosmic-ray acceleration and propagation in the Galaxy. 
Various interpretations have been put forward, such as 
the spectrum superposition of the local source and the 
back grou nd (T5| . the non-linear acceleration of cosmic 
rays and superposition of the injection spectra of 
the cosmic ray sources p^j . Instead of performing an ad- 
vanced study of the possible physical origin of the nuclei 
excesses, in this work we simply take such observational 
indication (i.e., the hardening of the cosmic ray proton 
and helium spectra above ~ 240 GV [l3|) as the main 
motivation to consider the possibility that the primary 
electron spectrum gets hardened at ~ 240 GeV. 

This work is structured as follows. In section II we take 
the latest AMS-02 positron ratio data to set new bounds 
on the dark matter models. In section III we discuss the 
possible positron-to-elcctron ratio signal in the presence 
of primary electron spectrum hardening and then com- 
pare with the observational data. The implication on 
the physical origin of the cosmic ray nuclei excesses is 
also investigated. We summarize our results with some 
discussion in Sec. IV. 
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II. AMS-02 POSITRON RATIO DATA: BOUNDS 
ON DARK MATTER MODELS 

Dark matter is a form of matter necessary to account 
for gravitational effects observed in very large scale struc- 
tures such as the flat rotation curves of galaxies and the 
gravitational lensing of light by galaxy clusters that can- 
not be accounted for by the amount of observed/normal 
matter [l^. The most widely discussed candidate is the 
so-called weakly interacting massive particles (WIMPs), 
which may annihilate with each other or decay and then 
produce particle pairs such as photons, electrons and 
positrons and so on p^ . That is why dark matter may 
be able to account for the observed positron and elec- 
tron excesses, as extensively examined in the literature 

[Hill. 

The high quality AMS-02 data confirms the positron 
excess discovered by PAMELA with unprecedented ac- 
curacy in the energy range of 10 — 350 GeV JJj. It is the 
time to explore the allowed WIMP parameter regions in 
both the annihilation model and the decay model. The 
electrons originated from dark matter annihilation/decay 
will suffer from inverse Compton scattering of interstel- 
lar background photons (e.g., cosmic microwave back- 
ground, dust emission and star light) and boost these 
photons to GeV energies, becoming part of the Galac- 
tic diffuse emission. Hence the Galactic diffuse emission 
detected by space telescopes can be used to reliably con- 
strain the physical parameters of dark matter particles 
p^ . The latest bounds set by the Fermi-LAT Galac- 
tic diffuse emission data are presented in [23 |. To make 
use of such bounds, our dark matter density distribution 
profiles as well as the cosmic ray diffusion parameters 
are taken to be the same as that adopted in [l^. The 
smooth dark matter density p are parameterized with a 
Navarro-Frenk- White (NEW) spatial profile 23] 

p{r) = Pod + r0/rs)V[(l + r/r,fr/rQ], (1) 

and the isothermal-sphere (ISO) profile ^] 

pir)^Poirl+r'J/{r'+rl), (2) 

respectively. Where po = 0.43 GeV cm^'^ is the local 
density of dark matter, tq = 8.5 kpc is the distance of 
the solar system from the Galactic center, ~ 20 kpc 
and Tc — 2.8 kpc are two scale radii. 

In the modeling, we consider the channels of annihila- 
tion/decay into iJi^ p~ and t+t^. These particles propa- 
gate in the Galaxy. In general, it is difficult to solve the 
propagation equation analytically, given the complicated 
distributions of the sources, interstellar matter, radiation 
field and magnetic field. In this work we adopt the GAL- 
PROP [1^ package to numerically calculate the prop- 
agation of the cosmic ray particles, including the con- 
tribution from dark matter annihilation or decay. The 
cosmic ray diffusion parameters are the following, includ- 
ing the Halo height — A kpc, the diffusion coefficient 
Dq = 5.3 X 10^^ cm^ s"\ the diffusion index (5=1/3, the 




FIG. 1: The regions of parameter space (blue and green re- 
gions, 99.5% and 99.999% confidence levels) which provide 
a reasonable fit to the positron fraction data of AMS-02 [111] 
comparing with the bounds set by Fermi-LAT Galactic diffuse 
emission (adopted from [22], the black lines represent the up- 
per limits on {av) found in the analysis with no model of the 
astrophysical background, while the red lines are the bounds 
found in the analysis with a modehng of the background). In 
each panel, the annihilation/decay channel as well as the dark 
matter distribution model have been marked. If one takes the 
bounds found in the analysis with a modeling of the back- 
ground of [S^l (i.e., the red lines), the channel XX ~^ A''''A'~ is 
allowed only for the isothermal-sphere like dark matter dis- 
tribution model and the decay channel x ^ M^M" is viable. 
We do not present the cases of annhilation/decay into t^t~ 
since they have been excluded (by the red lines). 

Alfven velocity Va = 33.5 km s~^, the nucleon injection 
indexes below and above the break rigidity pbr = 11.5 
GV are 1.88 and 2.39, respectively. 

In Fig|T] we present the parameter regions allowed by 
the latest AMS-02 data, where {av) (rdec) is the velocity- 
weighted-cross-section (lifetime) of dark matter annihila- 
tion (decay) into p^ pT . Our parameter space regions are 
much more compact than that obtained in f2lj since the 
AMS-02 data have much smaller errors and extend to 
much higher energies. Evidently, if one takes the bounds 
found in the analysis with a modeling of the background 
of 22] (i.e., the red lines), the channel xX ^ t^^l^-^ is 
allowed only for the isothermal-sphere like dark matter 
distribution model but the decay channel x ~^ M"*"/^" 
is still viable. The channels of annihilation/decay into 
T+r~ are not presented since they have been excluded 
by the Galactic diffuse emission data. 

Alternatively, the positron excess detected by 
PAMELA/Fermi-L AT/ AMS-02 may be mainly con- 
tributed by pulsars. It is widely assumed that the high 
energy electrons/positrons are generated through the 
cascade of electrons accelerated in the magnetosphere 
of pulsars [1^, [2^. The energy spectrum of e+e^ in- 
jected to the galaxy from pulsars can be parameter- 
ized as a broken power-law with the cutoff at Ec, i.e., 
dN/dE cx ApsrE'" exp{—E/Ec), where Ec ranges from 
several tens GeV to higher than TeV and the power-law 
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FIG. 2: The blue (red) solid line represents the best fit to 
the positron ratio data in the "background + pulsar" ( "back- 
ground + dark matter + pulsar") model. For the three- 
component model, we also show the positron-to-electron ratio 
(at energies > 240 GeV) expected in different primary elec- 
tron spectrum models: the solid line is for a smooth spectrum 
without hardening, while the dashed and dotted lines are for 
the spectra with hardening at the energy E — 240 GeV (the 
indexes of hardened primary electron spectra are marked in 
the Figure). The AMS-02 data are taken from 



index a ranges from 1 to 2.2 depending on the gamma-ray 
and radio observations (28l - [30| . Following 32] , E^. = 860 
GeV is adopted in the modelin g. T he spatial distribution 
of pulsars is parameterized as |31j 



/(i?,z)(x ( A 



exp 



b{R-RQ) 



exp 



(3) 

where Rq =8.5 kpc is the distance of solar system from 
the Galactic center, Zg « 0.2 kpc is the scale height of 
the pulsar distribution, a = 2.35 and b = 5.56. The 
other parameter appearing in the above equation is the 
normalization factor Apsr, and it will be determined 
in our modeling. Following (27j . we consider the two- 
component (background -I- pulsar) model as well as the 
three-component (background -I- dark matter + pulsar) 
model. The background electron injection indexes below 
and above the break rigidity pbr.c — 4 GV are taken to be 
1.54 and 2.6, respectively. For the two-component (three- 
component) model our best fit of the data at energies 
> 10 GeV gives a minimum = 53.35 (27.59), where 
the parameter a = 1.28 (1.265) is adopted. Please see the 
solid lines in Figl2]for the best fits. Such a result seems 
to favor the three-component model, in which the anni- 
hilation cross section < av >xx'^e+e-^ 10~^^ cm^ s"^ 
as well as the dark matter rest mass my, ~ 100 GeV are 
needed. As noted in , such a cross section is consistent 
with all of the current bounds of the indirect detection 
measurements. 



III. POSITRON-TO-ELECTRON RATIO 
BEHAVIOR: THE CASE OF 
PRIMARY-ELECTRON SPECTRUM 
HARDENING AT - 240 GEV 



ATIC, CREAM and PAMELA all found considerable 
cosmic-ray-nuclei spectrum hardening above the mag- 
netic rigidity ~ 240 GV. One would naturally speculate 
that the primary electron spectrum also gets hardened at 
~ 240 GeV. With the given electron/positron total spec- 
trum detected by (for example) Fermi-LAT, the harden- 
ing of primary electron component inevitably leads to the 
softening of the rest "additional component" which has 
been widely assumed to consist of electron and positron 
pairs. Hence the positron-to-electron ratio should be dif- 
ferent from that in the absence of a spectrum harden- 
ing. Such a possible modification is of our interest. The 
change of the spectral index of primary electrons at ^ 240 
GeV is hard to reliably estimate and we assume it is 
within the range 0.18 — 0.3, consistent with the proton 
(helium) spectrum hardening found by PAMELA p3| . 
The following calculation is straightforward. Firstly, for 
illustration, we take the analytical approach. As shown 
in the Fig. 23 of [265, the primary electron spectrum at 
240 GeV is oc E-^-^^ (see also FigH] of |271] for a rather 
similar result) and the total electron/positron spectrum 
is oc E'"^-^ for E > 240 GeV. Therefore, if the primary 
electron spectrum gets hardened by ^ ijo is (^e'^-^) at 
E > 240 GeV, the subsequent positron-to-electron ra- 
tio would get flattened (decreased)^. Secondly, to check 
our speculation more qualitatively, we perform a joint 
fit to the AMS-02 positron-to-electron data [33] and the 
Fermi-LAT electron/positron total spectrum data ^ 
and then evaluate the positron-to-electron ratio beyond 
the PAMELA energy range in the cases with and with- 
out primary electron spectrum hardening. The numerical 
results are presented in Figl^l As expected in our ana- 
lytical approach, the resulting positron-to-electron ratio 
at energies above ~ 240 GeV gets flattened or decreased, 
depending on the degree of the primary electron spec- 
trum hardening. The current AMS-02 data show tenta- 
tive evidence for flattening at energies above ^ 260 GeV. 
However the errors are still large and a robust conclusion 
is not possible. A full test of our prediction by AMS-02 
may be realized in a few years. On the one hand, the sim- 
ulation suggests that AMS-02 operating for 3 years will 
measure the positron-to-electron ratio at energy 400 
GeV with an error ^ 0.01 On the other hand, at 

400 GeV, the difference of positron-to-electron ratio in 
the cases with and without primary electron spectrum 



^ This conclusion is valid as long as the electron/positron total 
flux is similar to that measured by Fermi-LAT. If instead, the 
electron/positron total flux at energies 300 — 800 GeV resembles 
the ATIC data, it is diflicult to modify the electron-to-positron 
ratio unless the primary electron spectrum hardening is much 
harder than what we have assumed. 
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hardening can be as large as ^ 0.04 (see Figl2]). The 
detection prospect is thus promising. 

We would like to point out that the Australia Tele- 
scope National Facility (ATNF) pulsar model also pre- 
dicts fluctuations of the positron-to-electron ratio and 
sometimes the ratio drops off quickly (soj . However, un- 
less there is an additional hard electron component, the 
electron/positron total flux should display the same be- 
havior, which is different from our scenario. The same 
argument applies to the dark matter model as well. 

The spectrum hardening of both primary-electrons and 
nuclei at ~ 240 GV, if confirmed by AMS-02 in the 
future, is likely attributed to a "nearby" supernova- 
remnant-like source with a lifetime (in view of the sig- 
nificant softening of electron/positron total spectrum at 
the energy Ecoo\ ~ 1 TeV) 



10^^ s ( 



cool 



1 TeV 



Utot 



1 eV cm 



(4) 



wher 



MB+Ucmb + Mdust+Ustar, "B = /StT is the 



ere Mtot — "B "r"cmb"r "-dust ^ "-star, "B — ^IQ 

magnetic field energy density, Wcmb, "dust and Ustar are 
the photon energy densities of cosmic microwave back- 
ground, dust emission and the star emission, respectively. 
Correspondingly, the source should be within the char- 
acteristic cosmic ray diffusion radius 



1 kpc ( 



Dn 



E 



1028-7 cm^ s" 



where T) — Do{E/l GeV)* is the spatial diffusion coeffi- 
cient which is usually regarded as isotropic Q. In other 
words the source should be relatively nearby. There are 
two candidate supernova remnants, i.e., Geminga and 
Loop I 1131 . Here we do not consider nearby pul- 
sars since they are expected to produce electron/positron 
pairs rather than mainly electrons. The nearby super- 
nova remnants instead produce "ignorable" amount of 
positrons. This can be straightforwardly understood as 
the following. Let us make the "optimistic" estimate, i.e., 
assuming that the total cosmic ray protons above ~ 240 
GeV are dominated by the nearby source. The chance 
for one high energy proton to produce one positron is 
P ~ CTppncr/a ~ 3.5 x lO'^ (n/0.5 cm-3)(T/10" s), 
where a^p « 70 mb is the total cross section of production 
of pions in proton-proton collision, and n ~ 0.1 — 1 cm~^ 
is the number density of the local interstellar medium. 
The high energy proton loses about 20% energy in one 
proton-proton collision, roughly one quarter converts into 
positron via the decay of the positively charged pion 
(i.e., 7r+ /.t+ + 1^1^ e'^ + + Dfj, + v^). Hence 
at the energy of 240 GcV, the positron-to-proton ra- 
tio is - F/20i '^ - 2 X IQ-^ (n/0.5 cm-3)(r/10" s), 
well below the value ~ 3 x 10"'* inferred from the 
PAMELA/ AMS-02 data 0, [HI, where the E-'^ '^-Wke 
proton spectrum has been taken into account. If the 
future AMS-02 data show no evidence for the hard- 
ening of the primary electron spectrum at 240 GeV, 



the cosmic ray sources accounting for the nuclei ex- 
cesses likely locate at a distance further than R 
1.8 kpc ( " 



1028/ 



eV cm.""^ 



E 



24 TeV'' 



-1/3 



IV. CONCLUSION AND DISCUSSION 

The accurate positron-to-electron ratio data of AMS- 
02 confirms the positron excess discovered by PAMELA 
up to the energy ~ 350 GeV. The dark matter anni- 
hilation/decay can account for the positron excess but 
the Fermi-LAT Galactic diffuse emission data has im- 
posed a tight constraint on such a kind of possibility. 
The channels of annihilation/ decay into r+r^ pairs have 
been excluded. The channel XX ~^ M^A*^ is possible 
only for very cored dark matter distribution model (e.g., 
isothermal-sphere profile) while x — )• fJ-^fJ-^ is viable 
(see FigH]). Alternatively, the positron excess detected 
by PAMELA/Fermi-L AT/ AMS-02 may be dominated by 
pulsar-like astrophysical sources and the contribution by 
dark matter may be small but detectable. At energies 
< 100 GeV, a positron/electron component from either 
dark matter annihilation/decay or pulsar seems to be fa- 
vored by the AMS-02 data (see Figl2]). In the case of 
annihilation of dark matter with a rest mass ~ 100 GeV, 
the cross section into electron/positron pairs is needed 
to be ~ 10~2^ cm'^ s~^, consistent with all of the current 
bounds of the indirect detection measurements. 

Whether the primary cosmic ray electron spectrum 
gets hardened at ^ 240 GeV, like that observed in 
cosmic-ray-nuclei spectra, is an interesting question. If 
it does, the subsequent positron-to-electron ratio may 
get flattened or even decreased, depending on the degree 
of the primary electron spectrum hardening (see Figj^]). 
The current AMS-02 data show tentative evidence for 
flattening at energies above ~ 260 GeV. However the er- 
rors are still large and a robust conclusion is not possible. 
In the future, AMS-02 will precisely measure/separate 
the cosmic ray electrons and positrons up to ~ 1 TeV, 
with which the primary electron spectrum hardening hy- 
pothesis can be unambiguously tested. The spectrum 
hardening of both primary-electrons and nuclei at ^ 240 
GV, if confirmed by AMS-02 in the future, likely suggests 
a common physical origin. The electron spectrum soft- 
ening at ^ 1 TeV would thus point to "nearby" but old 
supernova-remnant-like source(s) with a lifetime < 10^'^ 
s, like Geminga and possibly also Loop I. 
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